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^sj We present results for the mass of the b quark in the MS scheme obtained by 

calculating the binding energy of the B meson in the static limit. The self energy 
of a static quark, Eq° needed for this purpose, is now known to 0(a s ) in the 
quenched approximation. We find a preliminary value of m^(m^) = 4.34(7) GeV 
at rtf = 0. The error is dominated by the remaining uncertainty in Eq°. In 
addition, using E^° at O(o 2 ), we estimate that the quark mass is reduced by 
approximately 70 MeV when two flavours of dynamical quarks are introduced. 
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1 Introduction 
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f— ^ . The mass of the b quark in the MS scheme (m&) can be extracted on the 

Q\ ' lattice using NRQCD, via the pole mass, M po i e , by calculating the binding 

energy of the B meson, Kund' M po i e — M^ p — Abi n d where M^ p is the 
spin-average of the experimental B and B* masses and Abind = E S i, n — Eq. 
E S i m is the energy of the B meson (at rest) in NRQCD and Eq is the b quark 
self energy. The pole mass is then converted to m^ at some scale /i using the 
continuum perturbative factor Z cont : m&(/x) = Z cont (^)M po i e . While M po i e 
O ' has an O(Aqcd) renormalon ambiguity, this is cancelled by similar effects in 

Zcont- mE is well defined. 

At present, Eo is only known to 0(a) for the b quark. However, in the limit 
of the b quark mass becoming infinite, Eq is known to 0(a 3 ) if internal quark 
loops are neglected (the quenched approximation). The tadpole- improved for- 
mula can be expressed as 



E^ = 1.070a p + 0.118a£ - 0.3(1.4)a$ : a p = 4™ /_0) (0.84/a). (1) 

The a p coefficient has been determined by Lepage et. al. 111. The error on 
the coefficient is numerical and quite large. However, it provides a realistic 
estimate of the uncertainty in Eq, compared to using 2 — loop perturbation 
theory and assuming the contribution of higher order terms is la p . 

In Eq. In Eq is expressed _in terms of a coupling constant defined on the 
lattice from the plaquette, a p a, evaluated at a characteristic gluon momenta, 
q* = 0.84/a, calculated using the BLM proceedureBu. In addition, the lattice 
calculation of E S i m has been tadpole improved, whereby all gauge fields on 
the lattice are divided by a 'mean-field' approximation to the gluon field, uq, 



to obtain more continuum-like operators. The corresponding tadpole improve- 
ment of Eq leads to the addition of the perturbative series for ln-uo- These 
ingredients result in a well behaved perturbative series for Eq. This is certainly 
not the case if the bare lattice coupling, oll = <7o/(47r) is used: 

££° = 2.1173a L + 11.152a! + 82.7(1.4)a| (2) 

Di Renzo et. al. U have also determined the 0(a 3 ) coefficient. They obtain 
86.2(.6), when at is used. Encouragingly, the two determinations, which have 
very different systematic errors, agree within 3<r. n 

Z cont has been calculated to <xt by Melnikov and van RitbergenQ: 



J cont 



1 - 0.4244a p - 0.4771c*p - 1.814a£ : a p = a£ 0) (0.62m£). (3) 



The series is well-behaved and we estimate the uncertainty in Z cont to be 3af,. 

The error introduced by working in the static limit, i.e. ignoring 0(Aqcd/M) 
contributions to Kundi leads to « 1% uncertainty in m&. The error arising from 
working in the quenched approximation is also likely to be around 1% (assum- 
ing a 10 — 20% shift in E S i m when sea quarks are included). These effects are 
the same size as the error arising from the numerical error in Eq. 

We obtained E^° m by extrapolating the simulation energy calculated at 
finite heavy quark mass. The latter was obtained as part of a high statistics 
study of the B meson spectrum in the quenched approximation at three lattice 
spaciiigs (a), with a -1 = 1—2.5 GeV. For details of the simulations see refer- 
enceLl Note that we use the spin-average of the experimental masses for the B 
and B* mesons in the expression for M po i e in order to reduce the error in using 
E^f m . In addition, we performed a study of sea quark effects using resu,Lf s ob- 
tained from a simulation including two flavours of sea quarks (n/ = 2) 
a -1 ~ 2 GeV. Only the O{oi) coefficient for n/ — 2 has been computed 
hence the comparison with m^ at ny = is performed using Eq and Z, 
this order. 

2 Results 

Tables |l| and |2j summarize our results. Within the combined statistical and 
systematic errors we see that this is the case and we take the result at j3 = 6.0 
as our best determination ofmi(mt) = 4.34(7) GeV at nj — 0. Note that the 
numerical error in Eq dominates the uncertainty in TrT£. 

In addition, using the results at /3 = 6.0 at nf — from the B s me- 
son and those obtained at nf = 2 we see that the b quark mass decreases 
by 70 MeV at O(oi) when sea quarks are introduced. Assuming the system- 
atic (perturbative) errors for the two simulations are correlated this is ~ 2a in 




Table 1: A Wnd and mtfmj) in GcV from the B meson at nj = 0. The statistical and 
main systematic errors are estimated, including those due to determining the inverse lattice 
spacing (a -1 ) and residual discretisation effects in E a i m ~ 0((\QCD a ) 2 )- 
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5.7 
6.0 
6.2 


.24 
.35 
.36 
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4.43 
4.34 
4.32 
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Table 2: The change in m b (m b ) from rif = to 2, where 2-loop perturbation theory is used. 



n/ mb(m6) stat. 



error 



4.45 
4.38 



.01 

.02 



the (remaining) statistical errors and the same size as the error in rnb(rfTb) at 
0(aZ). Further work is necessary to reduce the error in Eq, in order for sea 
quark and 0(Aqcd/M) effects to be significant. 
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